Tunneling conductivity in composites of attractive colloids 
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In conductor-insulator nanocomposites in which conducting fillers are dispersed in an insulating 
matrix the electrical connectedness is established by interparticle tunneling or hopping processes. 
These systems are intrinsically non-percolative and a coherent description of the functional depen- 
dence of the conductivity a on the filler properties, and in particular of the conductor-insulator 
transition, requires going beyond the usual continuum percolation approach by relaxing the con- 
straint of a fixed connectivity distance. In this article we consider dispersions of conducting spherical 
particles which are connected to all others by tunneling conductances and which are subjected to 
an effective attractive square well potential. We show that the conductor-insulator transition at 
low contents <j> of the conducting fillers does not determine the behavior of a at larger concentra- 
tions, in striking contrast to what is predicted by percolation theory. In particular, we find that 
at low </> the conductivity is governed almost entirely by the stickiness of the attraction, while at 
larger <j> values a depends mainly on the depth of the potential well. As a consequence, by varying 
the range and depth of the potential while keeping the stickiness fixed, composites with similar 
conductor-insulator transitions may display conductivity variations of several orders of magnitude 
at intermediate and large <j> values. By using a recently developed effective medium theory and 
the critical path approximation we explain this behavior in terms of dominant tunneling processes 
which involve interparticle distances spanning different regions of the square-well fluid structure as 
<j> is varied. Our predictions could be tested in experiments by changing the potential profile with 
different depletants in polymer nanocomposites. 

PACS numbers: 73.40. Gk, 82.70.Dd, 61.20.Ja, 64.60. ah 
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I. INTRODUCTION 

The challenge of understanding the electrical trans- 
port properties in conductor-insulator composites is cen- 
tral for conceiving and designing new composite mate- 
rials with unique electrical properties, and has fueled 
ongoing research ranging from fundamental statistical 
physics to the more applied materials science and nan- 
otechnology. In this respect, recent years have seen re- 
markable progress in the design and synthesis of poly- 
mer nanocomposites with controllable structural param- 
eters, so as to combine advantageous properties of the 
insulating polymer matrix (flexibility, light weight, trans- 
parency, etc.) with appropriate levels of electrical con- 
ductivity a. For example, high aspect-ratio conducting 
fillers, such as carbon fibers and nanotubes, graphite 
platelets and graphene sheets, as well as segregated dis- 
persions of carbon-black particles, have been shown to 
reduce strongly the loadings <j) needed to achieve con- 
ducting composites [3-0] ■ The interplay of multicom- 
ponent fillers has been recently demonstrated to have a 
strong influence on transport properties [4. Other stud- 
ies have evidenced the role of depletants Q and of shear 
forces Q in the formation and structure of the conduct- 
ing network of carbon nanotubes, thus underlining the 
importance of the coupling between the conducting and 



insulating phases Q, and how this can be used to tunc 
the overall composite conductivity. 

On the theoretical side, the general understanding of 
the effects of the filler properties on a has been tradition- 
ally based on the continuum percolation theory of trans- 
port. The premise behind this approach is that there 
exists a fixed inter-particle distance beyond which two 
particles are electrically disconnected, thereby leading to 
a well-defined percolation threshold <fi c for the connected 
particles which marks the transition between the con- 
ducting (at <j) > 4> c ) and the insulating (at 4> < 4>c) 
regimes of the composite [1,0. In the vicinity of 4> c , 
the resulting percolation conductivity is thus predicted 
to follow the power-law relation 
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where t ~ 2 is the universal transport exponent for three- 
dimensional percolating systems. Equation ([1} is custom- 
arily used to interpret the 0-induced conductor-insulator 
transition observed in real composite materials, which 
are thus viewed as truly percolating systems, although 
the values of t extracted from experiments often deviate 
from universality [l], [l(| . Due to the percolation hypoth- 
esis, most theoretical efforts have been concentrated on 
the calculation of the percolation threshold cj> c and of its 
dependencies on the filler particle shapes [TTl - fT3| , disper- 
sions |l^ , [l5l ]. an d filler- filler and filler-matrix interactions 

In polymer nanocomposites, or more generally in 
colloidal- like dispersions of conducting fillers in insulating 
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matrices, the electron transfer between particles is medi- 
ated by tunneling or hopping processes, which implies 
that the conductance between any two particles decays 
continuously with the inter-particle distance [see Eq. ^ 
below]. In this case, no fixed connectivity length can be 
unambiguously identified and, consequently, the premise 
for the existence of a percolation transition, and for the 
validity of Eq. ([T]), is unjustified for this class of com- 
posites. This poses the problem of understanding the 
observed filler dependencies of a without relying on the 
percolation hypothesis. 

Recently, the problem of finding the functional de- 
pendence of a for colloidal composites without imposing 
any a priori fixed connectivity distance has been tack- 
led by considering explicitly the electron tunneling con- 
duction between any two particles in the composite [20I — 
[22j . One important outcome was the realization that, 
despite the absence of a real percolation threshold, the 
transition from the conducting to the insulating regimes 
could still be characterized by a characteristic concentra- 
tion 4> p below which a is dominated by the (small but 
finite) conductivity a p of the insulating matrix. In par- 
ticular, by considering equilibrium dispersions of hard 
spheroidal particles, 6^ was shown to decrease with the 
particle aspect-ratio |20| , in accord with the general trend 
observed in composites with fibrous and plate-like fillers. 

In this article we study by theory and simulations the 
functional dependence of a when, in addition to steric 
interactions, the conducting colloidal particles are sub- 
jected to mutual attraction modeled by a square well po- 
tential. We show that attraction, when combined with 
tunneling, gives rise to unexpected features which could 
not be anticipated by the usual continuum percolation 
theory. Specifically, we find that the range and the depth 
of the attraction affect the conductivity in very different 
ways according to whether </> is small or large. In par- 
ticular, the low-0 behavior of er, and so the crossover (j> p 
to the insulating regime, turns out to be almost entirely 
governed by the overall stickiness of the attraction, while 



the conductivity at 



> 



is found to change by several 



orders of magnitude for different potential profiles even 
if the stickiness is kept constant. This kind of behavior is 
completely missed if the conductor-insulator crossover is 
treated as a true percolation transition occurring at the 
critical concentration (f> c , because in this case the 4> de- 
pendence of the percolation conductivity tx P erc of Eq. (TTJ) 
is basically determined by 4> c alone. 

As shown in the following, in which a recently devel- 
oped effective medium theory and the critical path ap- 
proximation are used to interpret the simulation results, 
the effects of the attractive potential on a can be ratio- 
nalized in terms of a ^-dependent dominant distance be- 
tween the fillers which governs the composite conductiv- 
ity. This correspondence establishes a direct connection 
between the functional dependence of a and the structure 
of the square well fluid, which could be in principle ex- 
ploited in the synthesis of real nanocomposite materials 
with optimized conduction characteristics. 



II. MODEL 

Attractive forces between filler particles in polymer 
nanocomposites may arise from several mechanisms such 
as van der Waals forces or depletion interactions induced 
by non-adsorbing polymers or surfactant micelles [23[ . In 
the latter case, the size and concentration of the dcple- 
tants control respectively the range and depth of the at- 
tractive interaction, which can be modeled by a suitable 
effective potential. Here, we represent the conducting 
fillers by hard-core spheres of diameter D and for any 
two spheres centered at and r^, we model the attrac- 
tion by a square- well potential u{rij) of the form: 



u(nj) 




for r%j < 
for D < ' 



D 



< XD 



> XD 



(2) 
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and XD (A > 1) is the range of 



attraction. In the following, we shall vary both A and 
e so as to consider the hard sphere (HS) case (e = 0), 
the short-range attraction regime e ^ and A < 1.25, 
and the adhesive hard sphere (AHS) limit [24[, which is 
obtained by taking A — > 1 and e — > oo with 



A exp(— e*) 
12(A- 1) 



(3) 



constant, where e* = e/kT and kT is the thermal energy, 
r is an inverse measure of particle stickiness which gov- 
erns, together with the concentration, the phase behavior 
of short-range square- well fluids [III [26| . In the follow- 
ing, we limit our study to values of r above the critical 
point for gas- liquid-like phase separation 27 1. 

We describe the electron transport processes by con- 
sidering each sphere as being electrically connected to all 
others through inter-particle tunneling conductances of 
the form 



g(rij) = soexp 



2(r; 



D) 



(4) 



where £ is the tunneling decay length, which is indepen- 
dent of the range of the attraction A. In writing Eq. (QJ 
we assume that the size D of the conducting particle and 
the temperature are large enough to neglect charging en- 
ergy effects and Coulomb interactions. Furthermore, we 
have ignored any dependence on r»j of the prefactor go, 
which we set equal to 1. 



III. NUMERICAL RESULTS FOR THE 
CONDUCTIVITY 

We performed standard Metropolis Monte Carlo sim- 
ulations to find equilibrium dispersions of N = 2000 
particles in a cubic box of side length L, which was 
changed to obtain different values of the volume fraction 
4> = ttND 3 /6L 3 . A dedicated Monte Carlo algorithm (28| 
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FIG. 1. (Color online) Conductivity probability function P(a) (thick solid lines) for various values of the volume fraction 
4> and for a square-well potential with A = 1.05 and r = 0.2. Here the tunneling factor is £/D = 0.01 and the number of 
conducting spheres is fixed at TV = 2000. The thin dashed lines are the Pcpa(o") curves obtained from P(S) of Fig. [5] and by 
using a = <to exp(— 25/£) with ao = 0.09. 



was used for the AHS limit. For all the 4> values consid- 
ered and for each potential profile we obtained Nr = 300 
independent equilibrium configurations of the system. 

To obtain the composite conductivity a, we con- 
structed for each of the Nr configurations a tunneling 
resistor network with inter-particle conductances g(rij) 
given by Eq. (UJ. In order to reduce computational times, 
we exploited the exponential decay of Eq. ^ by ne- 
glecting contributions from particles sufficiently far apart 
so as to reduce the number of connected particles in 
the network. To this end, we introduced an adjustable 
cut-off length S max such that the bond conductances for 
Sij > <5 max can be safely removed from the network with- 
out altering a [2(1 [H|. Finally, we solved this reduced 
network using a combination of numerical decimation 
and preconditioned conjugate gradient algorithms as de- 
scribed in Ref. I22I . From the resulting network conduc- 
tance G the dimensionlcss conductivity a follows from 
a = GD/L. 

In Fig. [1] we show some representative examples of the 
cumulative distribution function (CDF) P(a) of the con- 
ductivity (thick solid lines) obtained from all Nr real- 
izations of the network for A = 1.05, t = 0.2, and for 
different values of the volume fraction </>. In all cases the 
tunneling decay length has been set equal to£/.D = 0.01. 
Overall, the rise from to 1 of P(a) becomes more grad- 
ual as <p decreases, which comes from keeping the number 
of particles N fixed. Additionally, for <f> = 0.175 and 0.2 
the CDF displays a sudden increase at a ~ 0.5 x 10 ~ 7 . 
As discussed in more detail in Sec. [Vj this latter fea- 
ture can be traced back to the discontinuity of the radial 
distribution function of the square- well fluid at interpar- 
ticle distances r = XD. Since the numerical procedure 
for solving the tunneling network is rather time consum- 
ing, we have not attempted to study finite size-effects 



on P(<t) systematically. Instead, we have chosen to ex- 
tract the overall network conductivity from the condition 
P(a) = 1/2 applied to the CDF computed for N = 2000. 
As shown in Sec. [Vj where we study the conductivity 
within the critical path approximation for which a finite- 
size analysis is feasible, this criterion provides a robust 
estimate of the conductivity for infinite systems. 

In the main panels of Fig. [5] we show the network 
conductivity a (open symbols) as a function of for 
£/D = 0.01 and for different parameters A and e* of 
the square-well potential. For comparison, the HS limit 
is plotted in Fig. [2ja) as open squares. By inspection of 
the two panels of the figure we see that the </>dependcncc 
of er strongly depends on the way in which the potential 
profile is changed. Namely, in Fig. Uta) the conductivity 
is steadily enhanced for all </> values when the well depth 
is increased while keeping the range parameter fixed at 
A = 1.05, leading for example at <f> = 0.25 to an enhance- 
ment of up to 9 orders of magnitude compared to the HS 
limit. Conversely, as shown in Fig. HJb), when the po- 
tential depth is kept as fixed (e* = 1.0704) but its range 
varies from A = 1.05 to A = 1.25 the conductivity at 
low 4> is enhanced, while at intermediate and large <f> it 
is diminished by about three orders of magnitude. 

The different functional dependencies of a in Figs.^a) 
and (b) have important consequences in conjunction with 
the conductor-insulator transition of composites. Indeed, 
although we have considered conducting particles dis- 
persed in a perfectly insulating medium, in real polymer 
composites the conductivity a p of the polymeric matrix is 
small but finite at nonzero temperatures, being typically 
10 -13 — 10 -18 S/cm at room temperature. Hence, despite 
the fact that the tunneling conductivity drops all the way 
to zero as 4> — > 0, the conductivity of the complete sys- 
tem (particles and polymer) is limited from below by a p 
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FIG. 2. (Color online) Network conductivity a as a function of the volume fraction <f> obtained from MC calculations (open 
symbols) and from CPA (plus signs), (a) a for well range fixed at A = 1.05 and different well depths e* parametrized by r of 
Eq. ©. (b) cr for well depth fixed at e* = 1.0704 and A = 1.05, 1.15, and 1.25. Insets: corresponding EMA conductance g* . 



[29|. One can thus identify a crossover point separat- 
ing tunneling-dominated from polymer-dominated con- 
ductivity with the volume fraction <p p below which a 
matches that of the polymer matrix, as is generally done 
in the analysis of experimental conductivity data of real 
nanocomposites. As schematically shown in the main 
panels of Fig. [3J where the horizontal dashed line repre- 
sents a p fixed for illustrative purposes at 10 -20 , we can 
estimate <j) p as the intersection point between the calcu- 
lated tunneling conductivity and a p [2(ij . In this way we 
find that as the stickiness of the square- well is enhanced 
(i.e., t is diminished) <j) p is systematically lowered. How- 
ever, it is evident by comparing Fig. [^a) with Fig. [2£b), 
that the behavior of a in the conducting regime <f> > <fr p 
is not determined by <p p alone. For example the curves 
for r = 0.2 in Fig.^a) and for r = 0.22 in Fig.^b) have 
similar </> p but very different conductivities at larger <p. 

This latter feature is illustrated even more strikingly 
in Fig. [3] where a is plotted for different values of e* and 
A such that the stickiness parameter r remains constant. 
Even though <p p is basically unchanged in going from A = 
1.25 to the AHS limit at A = 1, the conductivity increases 
by 5 to 10 orders of magnitude above the HS case at fixed 
<j> in the whole <fi > <fi p range. Note that a similar behavior 
is also found if we change the value of a p provided that 
it remains sufficiently small. For example if in Fig. [3] we 
consider a p = 10 -15 , the curves for A < 1.15 would still 
have equal cf> p , though slightly larger than the previous 
case, but very different a at larger </>. 

Our general observation that in square-well fluids of 
conducting spheres the conductivity for <j> > <f> p is not 
determined by the position of the conductor-insulator 
crossover is in striking contrast with the usual contin- 
uum percolation description. Indeed if the conductor- 
insulator transition is viewed as a true percolation tran- 



MC 


1 

CPA 


i 






I 


— r- 


¥ 




— □— 

— o— 


+ HS 
+ li=1.2S 








^ 

v — 






- — A— 


+ A.=1 .15 


f 














— v— 


X=1.05 
















~ X = 


+ X=1.0 / 
0.2 J 


A' 

/ 
















i 


X 




10° 
















i<r s 
io-" 


















io-" 5 

10'° 


















1CT 5 
10™ 














I 




0.1 

1 


0.2 .0.3 
i 


0.4 0.5 



0.1 0.2 0.3 0.4 0.5 

<t> 

FIG. 3. (Color online) Network conductivity cr as a function 
of <j> obtained from MC calculations (open symbols) and from 
CPA (plus signs) for different well widths and depths with 
fixed t — 0.2. Inset: corresponding EMA conductance g* . 



sition, then the corresponding conductivity would follow 
Eq. |T]) for <f> > 4> c , where 4> c is the percolation threshold, 
and so the conductivity level in the conducting regime 
would be completely determined by (f> c . For the short- 
range potentials considered here 4> c is known to be re- 
duced by enhanced stickiness fl7| , in line with the behav- 
ior of c4 p , and consequently, from Eq. ([T]), the percolation 
conductivity for <fr > c4 c would also be systematically en- 
hanced. Conversely, different potentials with the same 
parameter r would give similar percolation conductivity 
curves, in strong disagreement with our results of Fig. [3] 
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FIG. 4. (Color online) Enhancement factor of the conductiv- 
ity for t — 0.2 compared to the conductivity in the HS limit 
for different values of the tunneling length £/D. 



The results shown in Figs. [5] and [3] also have important 
practical consequences in relation to the problem of con- 
trolling the conductivity of polymer nanocomposites by 
tuning the attraction through, for example, the depletion 
interaction Indeed, from Figs. [2] and 02 it turns out 
that the best strategy to obtain high levels of conduc- 
tivity for a broader range of filler concentrations is by 
choosing potential profiles that are very deep and short- 
ranged (but of course not strong enough to induce phase 
separation) . For depletion induced attraction, this corre- 
sponds to having a large concentration of depletants with 
small sizes compared to those of the conducting particles. 

We conclude this section by addressing the role of the 
tunneling decay length £ in the overall functional depen- 
dence of cr. In Fig. 2]we plot a for r = 0.2 and A = 1.05 
in units of the HS conductivity for £/D = 0.01, 0.02, and 
0.05. As expected, the strong enhancement obtained for 
(,/D = 0.01 is systematically reduced as £/D increases. 
For example at t;/D = 0.05 only a maximum 10 2 fold en- 
hancement with respect to the HS limit is achieved. This 
is because for large £/ D, two particles having separation 
less or greater than the potential range have similar tun- 
neling probabilities, and so the short-range structure of 
the square- well fluid is smeared out. If we consider that, 
depending on the material properties, the tunneling de- 
cay length ranges from a fraction of a nanometer to a few 
nanometers, from the results of Fig. [3] wc expect that par- 
ticle sizes of tens to hundreds of nanometers are optimal 
for a to be strongly enhanced by increased attraction. 



more complete understanding by employing the effective 
medium approximation (EMA) introduced in Refs. [2l| 
andUO- Briefly, this EMA amounts to replacing each con- 
ductance g{rij) of Eq. (|3]) by an effective value, g, which 
is independent of the inter-particle distance ry. By re- 
quiring that the resulting effective and fully-connected 
network has the same average resistance as the original, 
and by considering only two-site clusters 13011 . the follow- 
ing self-consistent equation is found [2l|, |30| : 



24<j) I dxx g 2 {x) 
'o 



1 



g* exp[273(x - l)/£] 



1 



2, (5) 



where x = r/D, g 2 (x) is the radial distribution function 
(RDF) for the conducting spheres, and g* — Ng/2 is the 
two-point conductance of the effective TV-node network. 

We have numerically solved Eq. for g* using the 
quasi-analytical model RDF proposed in Ref. [U for the 
square- well fluid, which reduces to the Percus-Yevick ap- 
proximation for both the HS and AHS limits and is in 
rather good agreement with the MC results of the RDF 
for A < 1.3 [U, [32 1 . Although g* is a two-point conduc- 
tance (which is the meaningful quantity for a complete 
network of identical resistors) rather than a conductiv- 
ity, the EMA results shown in the insets of Figs. [2] and [3] 
and in Fig. [5] (solid lines) closely follow the conductivity 
behavior obtained by the full numerical solution of the 
tunneling network. Using Eq. ([5]), we can thus explain 
the effect of attraction on <j by reasoning simply in terms 
of 172(2;) and of its dependence on the square- well poten- 
tial. To this end, let us rewrite the EMA conductance 
as 



9 



cxp[-2£(:r* - l)/£] , 



(6) 



where we have introduced a characteristic EMA distance 
x* . By noting that for £/D <C 1 the term in square 
brackets in Eq. (|5|) reduces to the step function 9(x* —x), 
x* is found as the solution of the following equation: 

24(f) [ dxx 2 g 2 {x) = 2, (7) 
Jo 

where the left hand side gives the number of particle cen- 
ters having distances less than x* from a particle at the 
origin. To proceed further, we note that for sufficiently 
narrow and deep potential wells the RDF can be approx- 
imated by 



g 2 {x)=e{x-l)[g 2 {l+)6(\-x) + l\, 



(8) 



where 172 (1 + ) is the RDF at contact, and from Eq. ((TJ) we 
obtain finally 



IV. EFFECTIVE MEDIUM THEORY 

The numerical results of Sec. IIHI and the effects of 
attraction on the functional dependence of a find a 



for x* < A, 



4^2(1+) 
-i--.g 2 (l+)(A 3 -l) + A= 



(9a) 



for x* > A. (9b) 




FIG. 5. (Color online) Comparison between the full numerical solution of the EM A equation (0 (solid lines) and g* obtained 
from Eqs. ((6]) and ((9]) (dashed lines). The parameters used in (a), (b), and (c) are the same of those in Figs. [2j a), Hfb), and [3] 
respectively. 



As shown in Fig. El g* obtained from Eqs. ((6]) and (|9]) 
(dashed lines) is in overall good agreement with the so- 
lutions of Eq. © with the full RDF of Ref. SH Further- 
more, Eq. © explicitly introduces two different regimes 
according to whether the EMA distance a;* is less or 
greater than A, which physically corresponds to whether 
the dominant tunneling processes occur between parti- 
cles with separation less or greater than the well range. 
The two regimes a;* > A and x* < A identify two different 
regions, respectively <j> < <fi* and (/)></>*, where <j>* is a 
characteristic concentration (not to be confused with <p p ) 
obtained from the condition x* = A applied to Eq. ([9|. 
i.e., 4(f)* g 2 (l + )(X 3 — 1) = 1. Since <?2(1 + ) scales essen- 
tially as exp(e*) [3l[, <fi* can be approximated by r/A 2 
for A — 1 and r small. 

From Eq. (|9aj) and (?2(1 + ) ~ cxp(e*) we thus see that 
the conduction at large densities {(j> > <j>*) is mainly gov- 
erned by the well depth, which is indeed expected since 
in this regime a;* < A. In contrast, from Eq. (|9b[) , g* 
turns out to be affected by both e* and A at low densities 
(</><(/>*) because the relevant tunneling distances span 
the entire well range (i.e., x* > A). In this regime and 
for sufficiently narrow wells, g2(l + ){\ 3 — 1) — A 3 reduces 
simply to l/4r — 1 + 0([X — l]/r), which explains the 
equal asymptotic behaviors for c/> — > of Fig. [3] as e* and 
A are changed but r is held fixed. 

From the low-density approximation derived above, 
it is also possible to deduce the conductor-insulator 
crossover point <j) p by introducing a two-point conduc- 
tance g p -C 1 of the polymeric matrix, which repre- 
sents the EMA equivalent of a p introduced in Sec. IIIII 
Thus by requiring that Eq. (|6]) coincides with g pi and 
using Eq. (|9b[) in the narrow well limit, we find for 
(£/2D)ln(l/ ffp )<l 



(10) 



Besides confirming that the structure of the square-well 
fluid affects <f> p only through the stickiness, Eq. (|10[) also 
makes it explicit that the locus of the conductor- insulator 
transition depends (although only logarithmically) on the 
insulating matrix through g p , in contrast to the usual 
continuum percolation theory which predicts a percola- 
tion threshold </> c independent of the matrix conductivity 

M. 



V. CRITICAL PATH APPROXIMATION 

Although colloidal dispersions of tunneling connected 
particles are intrinsically non-percolative systems, con- 
cepts and quantities of percolation theory are neverthe- 
less at the basis of the critical path approximation (CPA) 
[33| which, as already shown for HS fluids of conducting 
particles [2(| [22j , can reproduce to a high accuracy the 
tunneling conductivity behavior of composites. 

As we have seen in Sec. IIV1 for sufficiently small £/D 
the EMA conductance g* is dominated by a character- 
istic length, x* , such that the cumulative coordination 
number satisfies Eq. Q . In a similar way, CPA amounts 
to approximating the tunneling network conductivity by 



cr ~ Co exp 



26 c {</>) 



(ii) 



where <7o is a (^-independent conductivity prefactor and 
8 C (4>) = r c (<fi)~D is a critical distance given by the short- 



est among the Si 



D lengths such that the subnet- 



l + 6(£r/D)ln(l/ 5p ) 



work defined by all distances < S c (<p) forms a perco- 
lating cluster. From Eq. (fTTj) we see that CPA expresses 
a in terms of a critical connectivity distance, S c (4>), which 
is a genuine percolation quantity. However, contrary to 
the continuum percolation approach with fixed connec- 
tivity length, in CPA the critical connectivity changes as 
4> and the potential profile are varied. This adaptation 
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FIG. 6. (Color online) Percolation probability as a function of the connectivity distance 5 for various values of the volume 
fraction (f> and for a square- well potential with A = 1.05 and r = 0.2 (i.e. e* = 2.169). The number of conducting spheres is 
fixed at N = 2000. The vertical dashed line indicates 5/D = A — 1. 



of the connectivity range compensates for the artificial 
basis of the sharp cutoff at S c . 

To calculate 5 c ((f>) we follow the route described in 
Ref. HH. Namely, for fixed <fi an d potential profile, we 
coat each sphere with a concentric penetrable shell of 
thickness 5/2, and consider two spheres as connected if 
their penetrable shells overlap. A clustering algorithm 
described in Ref. [22| allows computation of the spanning 
probability P(S), which is plotted in Fig. [6]for A = 1.05, 
r = 0.2, N = 2000, and for different values of (j). Note 
that there is a rather sudden change in the slope of P(S) 
when 5/D crosses A — 1 (vertical dashed line) which is 
due to the discontinuity of the square well potential at 
r/D = A. 

In order to define a suitable criterion to extract 5 C from 
the P(S) curves for all parameter values used, and which 
is also insensitive to the feature at 5/D = A— 1, we 
have carried out the finite-size scaling analysis of the type 
shown in shown in Fig. where P(S) is compared with 
the percolation probability P{<j>) calculated as a function 
of for fixed 5. In Fig. 0a), we show the evolution of 
P(4>) as the system size L/D increases for A = 1.05 and 
5/D = A — 1 = 0.05. By using the finite-size scaling rela- 
tion 4> c — 4> C {L) oc L -1 /", where v ~ 0.88 is the correla- 
tion length exponent and <j> c {L) is such that P = 1/2, 
we find that the percolation threshold as L/D — > oo 
is C = 0.19035 ± 0.00001. In Fig. EJb) we plot P(5) 
obtained for different numbers N of particles by using 
4> = 0.19. Note that as N increases, the change of slope 
of P(S) at 5/D = A — 1 = 0.05 tends to decrease. This is 
more clearly seen in Fig. Etc), where we plot fits of P(8) 
with linear combinations of two error functions, and in 
the inset where we show the resulting probability distri- 
bution function (PDF) dP(8)/d5, which is characterized 
by two peaks of different strengths. If we use the cri- 
terion that 8 C (N) is given by P = 1/2, we find from 



5 C - 5 C (N) oc N~^ 3v that 5 C /D = 0.05066 ± 0.00007 
at N — > oo, practically coinciding with 5/D = 0.05 
of Fig. [7|a). Furthermore, for N = 2000 the criterion 
P = 1/2 gives 5 C / D ~ 0.0512, which is only 1% larger 
than the critical distance at N — > oo. 

From the finite-size study of Fig. [JJ and from other 
cases we have considered, we thus see that P(5 C ) = 1/2 
is a reliable criterion to find the critical distance even for 
N = 2000. When we compare in Fig. [1] the CDF of the 
network conductivity with Pcpa(c) obtained from P(5) 
by using the CPA expression a = (Tq exp(— 2<5/£) with 
(To = 0.09 (thin dashed lines), we see that the shapes 
of the two probability functions are very similar, which 
justifies the criterion P(cr) = 1/2 that we have adopted 
for the network conductivity results shown in Figs. @ 
and ©. 

Let us now discuss the conductivity dependence on <j) 
and on the potential profile in terms of the CPA formula 
(fTTT) . The calculated values of the critical distance are 
shown in Fig. [8] for the same parameter values as Figs. [2] 
and [21 For all cases studied, 5 C is a monotonically de- 
creasing function of </>, which is to be expected because 
larger particle concentrations bring about shorter mean 
inter-particle distances. However, the rate of decrease 
depends on whether 5 C /D is greater or less than A — 1 
(horizontal dashed lines in Fig. [5]), which reflects the ex- 
istence of the two different regions of the conductivity 
already discussed in Sec. IIVI In particular, as shown in 
Fig. [5(c), when 5 C /D > A — 1 the critical distance is gov- 
erned solely by the particle stickiness and closely follows 
the AHS limit, in full correspondence with the low-0 be- 
havior of a of Fig. [3] Indeed when the 5 C (4>) values of 
Fig. [S] arc inserted in Eq. ([TTj) . we find that the result- 
ing CPA conductivity (plus signs in FigsJ5] and [3]) is in 
excellent agreement with our MC values [34| . 

We therefore see that, in analogy with the results of 
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FIG. 7. (Color online) Percolation probability P for A = 1.05 and r = 0.2 (i.e. e* = 2.169). (a): P as a function of <j> for S/D 
fixed at A — 1 = 0.05 and for different sizes L of the cubic box. (b): P as function of S/D for 4> = 0.19 and for different number 
N of spheres, (c): fits of the percolation probability of (b) obtained by using a linear combination of two error functions. Inset: 
probability distribution function PDF obtained from dP(5)/dS. 
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FIG. 8. (Color online) Calculated critical distance S c as a function of (f> for (a): A = 1.05 and different well depths; (b): 
e* = 1.0704 and A = 1.05, 1.15, and 1.25; (c): A and e* varying with r = 0.2 fixed. The horizontal dashed lines identify 
5 C /D = A - 1. 



Sec. IIV1 the conductivity changes driven by attraction 
are due to the variations of the dominant tunneling dis- 
tances, which are represented by the effective connectiv- 
ity length S c . From Eq. (flT]) we see also that in the pres- 
ence of a finite conductivity a p of the insulating matrix, 
tunneling dominates the overall conductivity as long as 
<5c(0) < S c ((f> p ), where 5 c (<f> p ) = (£/2)/ ln(er /o>) is the 
connectivity at the insulator-conductor transition. 



VI. CONCLUSIONS 

In this article we have studied by theory and simu- 
lations the effects of short-range particle attractions on 



the functional dependencies of the conductivity a of equi- 
librium dispersions of conducting spherical fillers elec- 
trically connected through tunneling processes. Instead 
of employing the traditional continuum percolation ap- 
proach in which a fixed connectivity length is introduced 
to mimic the extent of the tunneling electron transfer, 
we have solved the network equations by treating each 
conducting particle as being connected to all others by 
tunneling conductances. We have shown that the lack of 
a fixed connectivity distance gives rise to new and impor- 
tant features which the continuum percolation approach 
fails to predict, and which could be crucial in the design 
of nanocomposite materials with unique electrical prop- 
erties. 
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In particular, we have found that the conductivity at 
intermediate and large contents of the conducting fillers is 
extremely sensitive to the profile of the attractive poten- 
tial, while the low-density regime is generally governed 
only by the stickiness of the attraction. As a conse- 
quence, materials with similar values of the conductor- 
insulator concentration <p p may display at <j> > <f) p con- 
ductivities which differ by several orders of magnitude if 
the attraction is tuned while keeping the stickiness con- 
stant. More generally, we have demonstrated that the 
knowledge of 4> p and of its dependence on the attraction 
are not sufficient conditions to presume the overall be- 
havior of a as a function of <fr, in striking contrast to the 
percolation transition approach represented by Eq. (fTJ) . 

As discussed in Sec. lIVI our simulation results can be re- 
produced by using a recently developed effective medium 
theory which, for sufficiently small values of £/D, ex- 
presses the conductivity in terms of the radial distribu- 
tion function at contact, explicitly relating the different 
regimes of conduction to the depth and range of the po- 
tential well. Furthermore we have shown how the critical 
path approximation overcomes the limitations of the con- 
tinuum percolation approach by allowing the connectiv- 
ity length to adjust for given <fi and potential profile, so 
to represent the dominant tunneling distances governing 
a. 

Although we are not aware of experiments on conduct- 
ing colloidal composites where the electronic conductivity 
is tuned by controlling the filler attraction [35| , our work 



may nevertheless stimulate experiments in this direction. 
In principle, our predictions may be systematically tested 
in experiments by changing the size and concentration of 
depletants in colloidal conductors. In particular we have 
predicted that an optimal enhancement of the conduc- 
tivity can be realized by high concentrations of deple- 
tants with small sizes compared to those of the conduct- 
ing particles, so to drive the system towards the AHS 
limit. Furthermore, stronger enhancements are expected 
in composites with conducting filler sizes sufficiently large 
compared to the tunneling decay length £, so to avoid a 
smearing out of the attraction potential well. 

We conclude by pointing out that the system consid- 
ered here, i.e., attracting fillers with spherical shapes, 
has served primarily to illustrate the general principles 
at the basis of the interplay between attraction and tun- 
neling. Stronger effects of the attraction arc expected 
when the conducting particles have shape anisotropy, as 
in fibrous or plate-like fillers, where the mutual orienta- 
tions of the particles couple with the attraction interac- 
tion. This represents certainly an issue of great interest 
for future investigations. 
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